Anomalous magnetotransport through reflection-symmetric artificial molecules 
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We calculate magnetotransport oscillations in current through a triple-quantum-dot molecule, 
accounting for higher harmonics (having flux period h/ne, with n an integer). For a reflection- 
symmetric triple quantum dot, we find that harmonics with n odd can dominate over those with n 
even. This is opposite to the behavior theoretically predicted due to 'dark-state' localization, but 
consistent with recent experiments [L. Gaudreau et al, Phys. Rev. B, 80, 075415 (2009)]. This 
feature arises from a more general result: In the weak-coupling limit, we find that the current is 
flux-independent for an arbitrary reflection-symmetric Aharonov-Bohm network. We further show 
that these effects are observable in nanoscale systems even in the presence of typical dephasing 
sources. 

PACS numbers: 73.63.Kv, 73.23.-b, 73.43.Qt 



The Aharonov-Bohm effect has been the subject of sus- 
tained experimental and theoretical investigation since 
it was first predicted in the middle of the 20 th cen- 
tury [TJ [S]. This effect is most commonly observed in 
transport through mesoscopic rings threaded by a mag- 
netic flux, in which the current or conductance is pe- 
riodically modulated as $ is varied [3H7]. For free elec- 
tron waves propagating in both arms of a ring, the pe- 
riodicity of these oscillations is typically the Dirac flux 
quantum, $o — h/e. However, weak localization due 
to disorder can also give rise to oscillations of period 
<i>o/2 [5] known as Altshuler-Aronov-Spivak (AAS) os- 
cillations. The interplay of these two effects has been 
studied theoretically [9] and controlled experimentally 
by embedding a quantum dot in one arm of the in- 
terferometer [TQ1 QT]. More recently, the fabrication of 
smaller and cleaner rings has allowed for the observation 
of higher harmonics — oscillations of period $0/^1 with 
n > 2 [HI [13] — associated with electrons circling the 
ring n times. Trajectories that circle the ring multiple 
times require longer coherence lengths to demonstrate 
robust interference and are therefore more susceptible to 
dephasing, typically leading to a decay of the harmon- 
ics with increasing n |12j . The smallest and simplest 
ring in which one could expect to observe these higher 
harmonics consists of a molecule of three sites in a trian- 
gular arrangement (provided, e.g., by quantum dots [14] . 
atoms, or implanted donor impurities |15j). 

Recent coherent magnetotransport measurements per- 
formed on a triple-quantum-dot device suggest, surpris- 
ingly, that the harmonics do not decay monotonically 
with increasing n [14]. Theory has predicted the forma- 
tion of a localized dark state in the triple dot, leading 
to a dominant n — 2 harmonic lfi 21 as in AAS oscil- 
lations 22 . In contrast, the results of Ref. [2] indicate 
that the n = 3 harmonic is dominant over that with 
n = 2. To the best of our knowledge, a detailed the- 
ory of the n > 2 harmonics in these systems has not 
yet been given. Because of their greater sensitivity to 
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FIG. 1. (Color online) Three quantum dots in a ring arrange- 
ment, with dots 1 and 3 connected to the source S and drain 
D via tunneling rates Fs and F_d, respectively. Dot i is sub- 
ject to a local potential Vi (taken here to be Vi = 0). The 
lowest single-particle orbitals of dots i and j are connected 
via tunnel couplings ry. The ring is threaded by a flux $. 
Inset: A reflection-symmetric Aharonov-Bohm network. 



electric and magnetic field fluctuations, understanding 
the large- n harmonics in small molecular rings may allow 
for enhanced magnetic-field and noise sensing in future 
nanoscale devices, similar to recent proposals for meso- 
scopic systems [53]. Moreover, in contrast with the case 
of mesoscopic rings, we show that discrete symmetries 
associated with finite-dimensional molecular systems can 
lead to non-trivial features in magnetotransport. 

For definitcness, we focus on the triple-dot arrange- 
ment shown in Fig. [I] although several results extend 
naturally to a larger number of sites and more general 
geometries. An accurate description of the large-n har- 
monics necessarily requires a theory that accounts for 
non- linear response (large voltage bias), that goes be- 
yond the leading order in weak-coupling to leads, and 
that accounts for strong Coulomb interactions. A co- 
herent master equation similar to that employed in Ref. 
[16] can be rigorously justified in an experimentally ac- 
cessible regime and satisfies all the above criteria. We 
show that for a reflection-symmetric configuration of the 
triple dot, the symmetry of the eigenstates forbids lo- 
calization that would lead to n — 2 AAS-like oscilla- 
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tions, consistent with recent puzzling experimental re- 
sults [14] . This effect arises from a more fundamental re- 
sult: In the weak-coupling limit (T/Ae —> with tunnel- 
ing rate T and molecular level spacing Ae), we find that 
the current through an arbitrary reflection-symmetric 
Aharonov-Bohm network (inset of Fig. [IJ is magnetic- 
field independent. 

In Fig. [T] we consider large single-dot level spac- 
ing, wo > Afi, U, U', with Afi = fig — Pd the bias 
for source(drain) at chemical potential Ps(d) and U, U' 
on-site and nearest-neighbor charging energies, respec- 
tively (setting e = H = 1). For simplicity, we take all 
dots to be at the same potential, Vi = 0, and choose 
Us > 0, pn < 0. In this regime only the lowest-energy 
dot orbitals (i = 1, 2, 3) participate. Orbitals i and j are 
connected via a tunnel coupling r^ , and dots 1 and 3 are 
coupled to the source and drain, respectively. 

We choose the gauge such that t\% and r 2 3 are real, 
while Ti 3 = |Ti 3 |e~ 27rl ^, with tf> = $/$o the reduced flux. 
The full Hamiltonian is 

n = n 3D + n L + n 3DL , (i) 

K 3D = T v d \a d 3<y + j Y( n it -nu.)+nc, (2) 

W.L = ^ e lkc\k<7 C lkcr\ Ti 3 DL = ^ Ulkd\ a Cl ka + Il.C. (3) 
Ika ilka 

Here, e z is the Zeeman splitting, the Coulomb interaction 
is He = UY Ji ni-\ni^+U'Y Ji> j^ <T in ie nj a ', d ia annihilates 
an electron on dot i with spin <r, and n Jcr = d\ a di a . The 
operator cika annihilates an electron in state k of energy 
eik and spin a in lead I = S, D. We neglect diamagnetic 
effects on \t^\, but account for the Peirls phase. 

We consider the Coulomb-blockade regime, ps < 
U, U' , in which energy-conserving transitions involve the 
vacuum and one-electron subspace. Restricting to one 
electron, the three localized orbitals hybridize into molec- 
ular orbitals with eigenenergies £/- CT for k — 1,2,3. We 
focus on the experimentally relevant high-bias regime, 
ps — Eka > k^T, Eka — Pd > &bT, where all one-electron 
states are accessible via sequential-tunneling processes 
with leads at temperature T. We further take the density 
of states, Vki, and dot-lead tunnel couplings, fyj., to be 
spin- and energy-independent over the bias Ap: vik ~ v\ 
and tuk ~ tn, leading to tunneling rates = 2Tri's\tis\ 2 
and T D = 2irv D \t 3D \ 2 . 

Dynamics giving rise to magnetotransport oscillations 
in this system are highly sensitive to sources of dephasing, 
so it is essential to establish an accurate regime of validity 
for the associated equation of motion, particularly in the 
chosen high-bias regime. To approach the problem sys- 
tematically, we start from the exact Nakajima-Zwanzig 
integro-diffcrential equation [24] for the full dot-lead den- 
sity matrix g(t) [IS]. From this, we find an equation for 
p, an effective reduced triple-dot density matrix in terms 
of spinless fermions. The dynamics of p are controlled by 



dot-lead correlation functions, which we evaluate within 
a Born approximation (leading order in Hzdl)- This ap- 
proximation is valid even for strong coupling, Ti ^> |ry|, 
as long as higher-order cotunneling processes can be ne- 
glected (see below). In the high-bias regime considered 
here, the dot-lead correlation time, r c ~ l/\pi — £fcer|, is 
taken to be much shorter than the characteristic evolu- 
tion time of p [> min (l/|ry |, 1/r;)]. In this singular- 
coupling limit [26, 27], a Markov approximation is justi- 
fied. We thus obtain a coherent master equation: 

p=-i [H, p] + 2T s V[d\]p + T D V[d 3 ]p. (4) 

Here, H = T ijd\dj is the triple-dot Hamilto- 

nian, where d\ creates a spinless fcrmion on dot i, 
d\ |0)=|i). The superoperator defined by T>[0]p = 
OpO^ - \ {0*0,p} is the Lindblad dissipator [21]. The 
first term in Eq. Q describes free evolution on the 
triple dot, while the two following terms describe tun- 
neling into and out of the triple dot. The factor 2 in 
the second term accounts for the two possible spin states 
into which an electron can tunnel. Contributions at 
higher order in W 3 dl can lead to dephasing and sup- 
press the coherent effects described by Eq. Q , especially 
in the considered high-bias regime. However, when < 
min dry |, Ys), we find that these cotunneling processes 
give a small correction in the considered regime, provided 

> max{(r s /A Ai ) 2 \ Tij \, [r^/tt/') 2 ] Ap}. 

The current is I = Tv{Ig}, where g = 
lim Tm ^ 00 (l/r m ) J Q m dtg(t) is the full stationary dot-lead 
density matrix with averaging (measurement) time r m . 
Here, the current operator is 2" = Md = i [H,Afo], where 
^Dkc c Dko- When cotunneling corrections are 
negligible, we find that I(4>) is given directly from the 
stationary population p 33 of dot 3, 

I(^=T D p 33 (^), I n = [ d^I^e 2 ™*. (5) 

Jo 

The behavior of the harmonics, /„, will generally de- 
pend on the choice of Tij. When all are equal, and 
for most generic choices of r^, the n-even harmonics, 
I^ki describing AAS-like oscillations, typically dominate 
in the absence of dephasing [TBT[22] . This behavior con- 
trasts with that suggested by recent experiments [14j . 
We have found, however, that the n-even harmonics can 
be suppressed for a range of tunnel couplings satisfying 
| rial — l T 23| 7^ | T-i 3 1 , which is indeed the case reported in 
Ref. [14]. As we will show, this suppression of AAS-like 
oscillations is a feature unique to low-symmetry molecu- 
lar systems and distinguishes this case from mesoscopic 
rings with a near-continuum of orbital states. To see this 
more explicitly, from Eqs. Q and (PH) we find a simple 
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analytical expression for / when t\i = T23 = r, |ri3| = t': 



I{4>) =V D 13 + z + 



[1 + x 2 ] [ysin(2^)+y 2 /2] 



1 + x 4 - 2x 2 cos(47T(/)) J 

x = t/t\ y = T D /r', z = T d /2T s . (6) 

The harmonics, /„, resulting from Eq. ([6]), are shown 
in Fig. [2](e) for typical experimentally realizable param- 
eters, displaying a non-monotonic dependence on n (in 
particular, \Is\ > I/2I), consistent with the experimen- 
tal findings [2H] of Ref. [HJ. Note that Eq. ^ does not 
satisfy the Onsager relation, i.e., I((f>) ^ I(— <j>), but the 
Onsager relation is generally not obeyed in the nonlinear 
(high-bias) regime considered here [2"§] . 

The harmonics /„ arise from paths that combine to 
encircle the ring n times [see Figs. |2ja)-(c)]. Trajecto- 
ries circling the ring multiple times become more sig- 
nificant when the dwell time on the ring ~ 1/Td is 
large compared to the timescale for coherent evolution 
1/Ae ~ 1/r', so it is useful to consider an expansion in 
y = T d /t'. The leading term, I ~ T D /{3 + z) oc y°, 
corresponds to the current from an incoherent (Pauli) 
master equation. The first subleading term (oc y) de- 
scribes harmonics with n odd and the next order (oc y 2 ) 
harmonics with n even. When Td/t' = y <C 1, the n-odd 
contributions dominate. The description of these higher 
harmonics requires that we go beyond the weak-coupling 
limit, i.e., beyond O(y ) [IJ. 

We note that for x = t/t' — > (resulting in ID trans- 
port through dots 1 and 3) the current, Eq. ([6]), still de- 
pends on cj) although there is only one path in this case, 
thus no flux enclosed. This seemingly unphysical result 
is a consequence of the non-commuting limits x — >■ and 
r m — > 00. We find that for any finite measurement time 
r m , I is ^-independent when x — > 0. Similarly, for any 
vanishingly small but finite x, the current / will have the 
(^-dependence indicated in Eq. ^ in the limit r m — > 00. 

In general, we find it convenient to describe magneto- 
transport oscillations in terms of two distinct contribu- 
tions. The first, static contribution, gives rise to AAS- 
like oscillations due to localization and derealization 
of molecular eigenstates as the magnetic field is varied. 
These static contributions can be found accurately from 
an incoherent (Pauli) master equation, which neglects 
off-diagonal elements of p in the eigenbasis of H and 
gives magnetotransport oscillations at leading order in a 
conventional weak-coupling expansion, O(y ). The sec- 
ond, dynamical, contributions give rise to the standard 
Aharonov-Bohm effect and its harmonics associated with 
the coherent motion of electrons circling the ring. The 
dynamical contributions arise only at subleading order in 
a weak-coupling expansion [~ 0(y) or higher]. 

We can understand the unusual behavior shown in 
Fig. [2je) by accounting for both the static and dynami- 
cal contributions. The static contributions are best ana- 
lyzed within degenerate perturbation theory |22) . When 




FIG. 2. (Color online) (a)-(c) Trajectories leading to the 
n — 1,2,3 harmonics. For n even, the transit times for the 
two paths are always different, suppressing the associated in- 
terference, (d) Eigenenergies e of the triple dot as a function 
of the reduced flux (f> for r = T12 = T23 = 1 7-13 1 (dashed) 
and for r = T12 = T23 7^ |t"i3 | = r' (solid). A gap opens at 
values of <j> — m/2 with m £ Z, protecting the state from 
localization, (e) Magnitude |j n | for r = 30 ^ieV, r' = 20 fJ,eV, 
2Ts = 20/ieV, and Td = 1 /LteV. The n-even harmonics are 
strongly suppressed with respect to those with n odd. 



the couplings have equal magnitude, |ryj = r the 
eigenstates are completely delocalized molecular states, 
pairs of which are degenerate for <j) = m/2, m € Z 
[see Fig. ^d)] 18J. Introducing a small real perturba- 
tion to the triple-dot Hamiltonian taking it away from 
= r Vz,j, leads to eigenstates that are symmet- 
ric and antisymmetric linear combinations of the de- 
generate states at leading order in degenerate perturba- 
tion theory. These combinations are typically localized, 
strongly suppressing current at near-degeneracies, when 
4> = m/2, and hence inducing AAS-like oscillations of 
period $ /2 [I6H201 [22] . However, when r = t 12 = t 23 
and t £ t' = \t u \, i.e. when [H, 9] = 0, 9 = UK 
being an anti-unitary operator composed of the parity 
operator II = |1)(3| + |3)(1| + |2)(2| and of the complex- 
conjugation operator K, we exploit this special discrete 
symmetry of the triple dot [30] to gain further insight. 
In this case the eigenstates \e p ) = J2i c p I*) 01 H are s i~ 
multaneous eigenstates of 9, uniformly delocalized across 
dots 1 and 3 (|cj| = |c||Vp), strongly suppressing the 
AAS-like oscillations. Indeed, a state |e p ) is loaded at 
a rate 2j$ = 2|c p | 2 rs and unloaded at a rate j D = 
|c„| 2 r£>. Solving the steady-state Pauli master equation 



2jgp 00 — 7^>Ppp = to obtain the 'static' contribution and 
using \cp\ 2 — \Cp\ 2 , we find that p pp /poo — 2T S /T D is the 
same and independent of (f> Vp. Because J2 P Ppp = ^ 
this means that the current, I — ToP33, is also inde- 
pendent of (f>. We emphasize that this ^-independence 
in the weak-coupling limit, T^/Ae — > 0, holds for any 
Aharonov-Bohm network with a single dot connected 
to each lead whenever the dot configuration has mirror 
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symmetry (inset of Fig. [I]). Thus, when [H, 0] = 0, 
any Aharonov-Bohm oscillations must arise from coher- 
ent dynamics on the triple dot contained in the first term 
on the right-hand side of Eq. Q and not from the local- 
ization/delocalization of eigenstates responsible for AAS- 
like oscillations. If a perturbation V breaks parity at 
4> = m/2 ([V,II] 7^ 0), lowest-order perturbation the- 
ory shows that the eigenstates will remain delocalized as 
long as | (e + |V|e_) |/A <C 1, where |e + ) and |e_) are the 
nearly degenerate eigenstates of H and II at <j> = m/2 and 

where A = 3|r'| — \Z|t'| 2 + 8r 2 /2 is the gap opened at 

the degeneracy point [see Fig.[2|d)] . Even when the AAS- 
like oscillations are suppressed, we might still expect to 
find a significant dynamical contribution to the n = 2 
harmonics coming from the interference of paths circling 
the ring twice [see Fig.J^b)]. However, because backseat - 
tering is forbidden in the high-bias regime, the remaining 
asymmetric paths result in different transit times, reduc- 
ing interference. In contrast, the n-odd contributions can 
arise from paths that are symmetric in both arms [see 
Figs. |2j>) and (c)]. 

To account for dephasing we add a term Vg(<) = 
— e^ i(T E(i) • Vin i(T to H, Eq. 0. This term describes 
electric-dipole coupling of localized dot orbitals at posi- 
tions to a fluctuating electric field, E(t), taken to be 
uniform across the triple dot. The effect of Ve(£) is to 
introduce a dissipator to Eq. Q: 



P . 



(7) 



i>3 



where Co generates the right-hand side of Eq. Q . When 
Vg;(i) dominates over % and when the noise is Gaussian 
and white, we find Tf^ ~ [r, — rj | 2 , implying that electric- 
field-induced dephasing is more significant for larger sys- 
tems. Note that Eq. Q is only strictly valid when Tf^ = 
or ^> max(ry,r;). When the dephasing term in 
Eq. ufh is non-zero, it acts as a which-path measurement 
for both arms of the ring and destroys the interference 
that would lead to higher harmonics. Fig. [i] shows \I n \ 
for the simple case where Tfj =Te for all Increased 

dephasing re-establishes the monotonic decrease of |/ n | 
with n when 3> i.e., when the characteristic 

evolution time | j | 1 on the triple dot is longer than the 
dephasing time r^ -1 . Thus, we do not expect the non- 
monotonic behavior of |/„| displayed in Fig. [ija) to be 
observed in larger rings where the relative electric dipole 
moment ~ |r< — Tj\ between arms is much larger. 

In summary, we have derived a non-linear (high-bias) 
transport theory for a triple quantum dot in the strongly- 
interacting Coulomb-blockade regime, and have used it to 
describe $o/n harmonics of Aharonov-Bohm oscillations 
beyond the limit of weak coupling to the leads. We have 
shown that for a reflection-symmetric Aharonov-Bohm 
network, the current is flux-independent in the weak- 
coupling limit (r/Ae —> 0). In the simplest case of a 
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FIG. 3. (Color online) Magnitude |J n | of the n = 1,2,3 har- 
monics as a function of the dephasing rate, Te, for Tvz = 
r 23 = 30/ieV, |n 3 | = 20 peV, 2r s = 20^ieV and Y D = lfieV. 
When r_E 3> \Tij\, \I n \ decreases monotonically with n. The 
first vertical line (at r_E = 1GHz), corresponds to the mea- 
sured dephasing rate for a nanoscale double quantum dot of 
size ~ 250 nm [31]. The second vertical line (at Fe = 16 GHz) 
corresponds to the estimated dephasing time for a meso- 
scopic ring of size ~ 1 fim |12| and is obtained from our 
theory by scaling according to the relative dipole moments 
[16GHz = (l/im/250nm) 2 x 1GHz]. 



reflection-symmetric triple dot, this results in a strong 
suppression of the n-even harmonics compared to those 
with n odd. This happens for two reasons: first, the 
n-even AAS-like harmonics due to the localization of 
triple-dot eigenstates are suppressed since symmetry re- 
quires that the eigenstates be delocalized across the triple 
dot; second, the dynamical paths leading to n-even har- 
monics have different transit times, suppressing interfer- 
ence. These results may help explain the experiments 
of Ref. [14] . We have also shown that strong dephasing 
restores a monotonic suppression of harmonics with in- 
creasing n. Finally, we remark that these results are im- 
portant for the design of noise and magnetic-field sensors 
using Aharonov-Bohm interference [53] . The coherent ef- 
fects presented here are likely to be substantially stronger 
in still smaller devices, such as those very recently con- 
structed from idividual atoms or donor impurities 15 . 
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